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ABSTRACT 

In the first of a series of papers investigating emission from blazar jets from radio to high- 
energy gamma-rays, we revisit the class of models where the jet has a uniform conical ballistic 
structure. We argue that by using simple developments of these models, in the context of 
new multi-frequency data extending to 7-ray energies, valuable insights may be obtained into 
the properties that fully realistic models must ultimately have. In this paper we consider the 
synchrotron and synchrotron-self-Compton emission from the jet, modelling the recent simul- 
taneous multi -wavelength observations of BL Lac. This is the first time these components have 
been fitted simultaneously for a blazar using a conical jet model. 

In the model we evolve the electron population dynamically along the jet taking into 
account the synchrotron and inverse-Compton losses. The inverse-Compton emission is calcu- 
lated using the Klein-Nishina cross section and a relativistic transformation into the jet frame, 
and we explicitly show the seed photon population. We integrate synchrotron opacity along the 
line of sight through the jet plasma, taking into account the emission and opacity of each sec- 
tion of the jet. In agreement with previous studies of radio emission, we find that a conical jet 
model which conserves magnetic energy produces the characteristic blazar flat radio spectrum, 
however, we do not require any fine-tuning of the model to achieve this. Of particular note, in 
our model fit to BL Lac — which at ^ lO'^^Ty is a relatively low jet-power source — we find no 
requirement for significant re-acceleration within the jet to explain the observed spectrum. 

Key words: Galaxies: jets, galaxies: active, radiation mechanisms: non-thermal, BL Lacertae 
objects: individual (BL Lacertae), radio continuum: galaxies, gamma-rays: galaxies. 
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1 INTRODUCTION 

Blazars are the most luminous and highly variable of active galac- 
tic nuclei (AGN). They have been observed at all wavelengths from 
radio to high energy 7-rays and are believed to be AGN observed 
at small angles to the jet axis ( jUrry & Padovani] (fl995}) resulting 
in Doppler boosting of their emission. The characteristic spectral 
shape is a non-thermal continuum with two main peaked compo- 
nents; the radio through to UV/x-rays show signs of polarisation 
typical of synchrotron radiation, the second component extends 
from low energy x-rays through to very high energy 7-rays and is 
normally attributed to inverse-Compton scattering of photons from 
the synchrotron-emitting electrons. The radio emission of blazars is 
predominantly observed to be flat or nearly flat in flux density (see 
for example [Owenet al.| ( 1980) and jAbdo et al.| ( |2010^ ) whilst the 
majority of the synchrotron and inverse-Compton emission can be 
well described by a power law and a turnover. 

In general, VLBI images of AGN observed at large angles to 
the jet axis show a continuous axisymmetric plasma jet, which is 
approximately conical in geometry, (see e.g.,|Kovalev et al!|(|2007|) 
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and [Krichbaum et al.| J2006l) with evidence for a small degree of 
curvature close to the base of the jet jAsada & Nakam ura (20lT|) 
and a blunt base with a wide opening angle ( |Hada et al.| ( |201l| ). 
Recent work by [Sokolovsky et al.| pOTTTl measured the frequency 
dependent core-shift for 20 AGN and found the results favoured a 
conical jet which conserved magnetic energy. 

The early and influential model of |Blandford & K6nigl| ( |1979] ) 
used a uniform conical structure for the jet, and was able to pro- 
duce the flat radio spectrum if magnetic energy in the jet was con- 
served. Further investigations by |Marsche r ( 1980), Konigl ( |1981^ , 
Reynolds ( 1982), M ufson et al . 1 1984) and Ghi sellini et al.| ( |1985| ) 
indicated that this class of model was a good fit to the observed syn- 
chrotron emission. However, at the time, 7-ray observations were 
rare, and replenishment of adiabatic and radiative losses in some of 
the models were deemed to be problematic. 

More recently, models for the radio and inverse-Compton 
emission from relativistic jets have predominantly considered all 
the emission to come from a small region of relativistic plasma, and 
such models have the ability to match the 7-ray spectrum of flares 
in blazars, e.g.|Kirk et al.|([T998l),|B6ttcher & ChiangH2002) , p^ 



|Kusunose| ( |2000 1 and Tsang & Kirk i 2007^ . The spherical blob is 
injected with a power-law (or broken power-law) electron energy 
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spectrum at regular time intervals and the distribution of electron 
energies is evolved using the continuity or Fokker-Planck equa- 
tion and the resulting emission spectrum is averaged over time. The 
continuity/Fokker-Planck equation takes into account the change in 
the electron population due to electron energy losses from emis- 
sion and due to the injection of additional electrons, but the energy 
injection mechanism remains arbitrary. 

Models which contain some extended jet component explicitly 
include those of, e.g.,|Marscher & Gear|(|1985|), Kaiser et al. (2000^, 
[Markoff et aT (20001, Spada et al. (2001 1 and Jamil et al. (2qTo}. 
These have variously included some attempts to treat particle re- 
acceleration, but are often aimed at x-ray binary systems and have 
not been concerned with the entire spectrum of a blazar from radio 
to 7-ray . Recently the conical jet was investigated by |Kaiser| ( |2006[ ) 
as part of an analytic study of the range of conditions in the jet 
plasma which could result in a flat radio spectrum, although this 
model too did not investigate the inverse-Compton component of 
the spectrum. 

With the advent of high-quality simultaneous radio-to-7- 
ray spectra, especially including data from the Fermi satellite, we 
were motivated to re- visit the problem of the full blazar spectrum. In 
particular, we wish to stress that the ultimate goal should be a model 
which contains both explicit and realistic jet parameters along with 
a detailed treatment of the particle emission processes and energet- 
ics within the jet. 

In this paper we set out the details of a flexible and rigorous 
model which will form the basis of a series of investigations into 
the properties of blazars. We use a simple conical jet geometry as 
a first-order model for the structure of the jet. Our focus is on ac- 
curately calculating both the synchrotron and the synchrotron-self- 
Compton components of the jet spectrum. We do this using a full 
Klein-Nishina calculation of the inverse-Compton scattering and in- 
tegrating the synchrotron emission along the line of sight through 
the relativistic jet, and we compare with the recent observations of 
BL Lac from radio to 7-ray . For a conical jet this is the first time a 
full-spectrum fit to a blazar has been attempted, and we also recover 
the jet opening angle, equipartition fraction and bulk kinetic power 

In this first paper we investigate whether a simple ballistic con- 
ical jet can reproduce the spectrum of BL Lac without need for 
significant acceleration or external Compton radiation. We include 
radiative electron losses but we do not yet include adiabatic losses 
in the jet. For the relatively low-power sources such as this, we find 
that there is not a need for significant re-acceleration within the jet. 
In the next paper in the series we analyse whether this is the case 
for Compton-dominant high energy sources. 

Most previous investigations of the conical jet have used sim- 
plifying approximations when calculating the opacity and observed 
synchrotron emission at small angles through the jet and approxi- 
mate analytical expressions for the inverse-Compton emission. We 
treat the calculations of the line of sight synchrotron opacity and 
inverse-Compton emission thoroughly. We use numerical integra- 
tion of exact expressions to calculate the inverse-Compton emis- 
sion. Unlike all previous investigations we also take into account 
electron energy losses along the jet due to both synchrotron and 
inverse-Compton emission and take into account the effect of these 
energy losses on the line of sight synchrotron opacity. 

In this paper we will first introduce the parameters and as- 
sumptions of our jet model including its geometry and electron pop- 
ulation. We then calculate the inverse-Compton spectrum produced 
from the scattering of an arbitrary isotropic seed photon population, 
although we shall only consider SSC seed photons in this work. We 
calculate the observed synchrotron emission from the jet and the 



synchrotron photon population. We show that the ballistic conical 
jet model is consistent with recent multi-frequency observations of 
BL Lacertae. We show that this model reproduces the observed flat 
radio spectrum if the magnetic energy is conserved along the jet 
and that this flat radio spectrum is a property which is insensitive to 
the model's parameters. Finally, we explicitly show the seed photon 
distribution and the evolution of the electron energy distribution. 



2 JET MODEL 

As described above, VLBA images of AGN indicate that jets flows 
are axisymmetric with a constant opening angle over much of the 
jet length and so can be modelled by a conical geometry. This ge- 
ometry was also favoured by |Sokolovsky et al.| ( |20Il"| ) to explain 
the frequency dependent core-shifts they measured. We will assume 
that the jet has an opening angle ^opening, initial radius Rq and 
length L, all defined in the centre of momentum frame of the fluid 
(fluid frame), see Figure IT] We assume that the jet is composed of 
an electron-positron plasma (for simplicity we will use 'electrons' 
to refer to both the electrons and positrons) which is injected at 
the base of the jet and the electrons lose energy radiatively as they 
move along the jet. Further, we make the assumption that both the 
electrons and magnetic field are isotropically and homogeneously 
distributed within the fluid rest frame whilst the fluid is moving in 
the lab frame with an associated bulk Lorentz factor 7buik. 

We assume that the jet structure is constant with time in the 
lab frame and that the geometry of the jet is that of a truncated 
cone of length 7buiki where the length of the jet in the lab frame 
is related to that in the fluid frame by a simple Lorentz contraction. 
We define the variable x as the length along the jet axis (in the fluid 
frame), where x = is defined dynamically as the base of the jet 
and X = L a& the end of the jet in the fluid frame. The definition 
of X is dynamic because in the fluid frame the structure of the jet 
is moving whilst the plasma is stationary, as opposed to a moving 
plasma and stationary jet structure in the lab frame (simply due to 
the relative velocity of the plasma and jet structure). We can now 
parameterise all of our physical quantities as functions of x. We 
also define the radius of the jet as R{x) where 



R(x) = Ro + a;tan6loponmg. (1) 



2.1 Determining the Jet Parameters at the Base 

We consider a jet with total power in the lab frame Wj . At the 
base of the jet the jet material has speed /3c (where we assume a 
relativistic flow P ~ 1). So the energy contained in a section of 
plasma of width one metre in the x-direction in the lab frame is Ej 
which is equal to — -. If the jet has a bulk Lorentz factor 7buik then 
the energy Ej in the lab frame is related to the energy contained in 
a one metre section of plasma in the frame of the fluid £", ""* by 



E^ 



1 ' 
7bulk 



E, 

7bulk 



t/e + [/fl, 



(2) 



where f/e and Ub are the electron energy and magnetic field 
energy contained in the first metre of the jet in the a;-direction, in 
the fluid frame. We define the equipartition fraction Aoqui as 



.'T-cqui — 



Vb_ 



(3) 
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Figure 1. Schematic diagram of the conical jet, showing parameters defined 
in the fluid centre of momentum frame 



function of the radius of the jet so that the total magnetic energy is 
conserved in a segment. 



Ub = nR{x 



,2B{xf 



B(x) 



2fioUB 



2^0 ' ' ' V TvR{x)^ • 

Substituting in the value of Ub from equation|9]we find. 



(9) 



B{x) = Bo 



Ro 
R{x)' 



(10) 



In this case the local magnetic energy is conserved as well as 
the radial and toroidal magnetic flux, however, the magnetic flux 
parallel to the jet axis increases with x. If the magnetic field is ran- 
domly oriented on small scales then the total flux parallel to the jet 
axis averages to zero, so the magnetic flux parallel to the jet axis 
does not change along the jet. 



Let us now consider the electron population in more detail. It is 
believed from optically thin synchrotron observations of AGN that 
the electron energy distribution is described by a single power law. 
We impose an exponential cutoff after -Emax i.e. 



N4E, 



ae: 



(4) 



where N^{Ee) is the initial injected electron distribution con- 
tained in the first metre of the jet, A and a are constants and a is 
thought to be between 1 and 3 from the theory of shock acceleration 
( |Beiil ( [T978l l, |Bell et al.| ( |20TT) and |Summerlin & Baring] pOTT] l). 
We assume that the electron population within a section of length 
Ax is independent of x. 

We define Bq to be the magnetic field strength at the base of 
the jet in the lab frame, at the base of the jet Ue and Ub are given 
by 



J E, 



U, = I E,N4E,)dE, 



Ub=7tRI^, (5) 



Ue = 



Ej 

7bulk(l + ^cqui)' 



^s = -^ T-. rrr^- (6) 



We can use these equations to calculate the coefficient A and 
to find a relation between Bo and Ro. 



A 



(2-Q)t/e 



(2 - a)E, 



(E, 



2-cv 
max 



Ei-r^) 7Llk(l + Aqni)(i5i 



E^:)' 



Ro 



^UBfJ-O 
■KlcBl 



2_B_,ylequiM0 



(7) 



(8) 



We choose Bo to be our independent parameter and deter 
mine Ro from the above equation. For reasonable values of Bq < 
10 ~^T, Wj > IQ^^ W and 7buik < 100 ( JTavecchio et al.| 
and|Lobanov](|2010|>), we find i?o > 3 x lO^^m. 
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2.2 Magnetic Field Strength 

For a jet with a constant bulk Lorentz factor which conserves mag- 
netic energy in each segment the magnetic field will change as a 



3 ELECTRON RADIATIVE LOSSES 

In our calculations we evolve the electron population dynamically 
along the jet by taking into account energy losses from synchrotron 
and inverse-Compton radiation. For simplicity we will consider the 
electrons in a slab of width one metre in the fluid frame as it propa- 
gates along the jet (or more correctly, since the slab is at rest in the 
fluid frame, the jet structure moves relative to the slab). As a section 
of the jet of width dx in the fluid frame moves past an expanding 
stationary slab, the electrons in the slab loose an amount of energy 
equal to the total power emitted by the entire section in the fluid 
frame divided by c, independently of the width da; of the section. 
This is because at any point the section contains n = dx one metre 
slabs and the end of the section (which is moving in the fluid frame) 
travels with speed c towards the slabs, so a slab takes " seconds to 
cross the section in the fluid frame. The evolution of the electron 
population along the jet due to energy losses can then be calculated 
using the equation below. 



iVe(-Ee,a; + da;) = Ne{Ee,x) 



Ptot{x,dx,Ee 



(11) 



where Ptot is the total power emitted by electrons of energy 
Ee by a section of jet of width da; in the fluid frame due to syn- 
chrotron and inverse-Compton processes. This equation is equiva- 
lent to solving the Fokker-Planck equation for an electron distribu- 
tion contained in a section of plasma which loses energy radiatively 
whilst moving along the jet. We can take into account the dynamic 
electron population by modifying the variable A to include a de- 
pendence on jet length and electron energy A{x, E^), which takes 
into account the changes to the electron population as a function of 
jet length and electron energy. 



A{x,Ee) = Ax 



Ne{Ee,x) 

N,{Ee,X = 0)' 



(12) 



We can now calculate the synchrotron and inverse-Compton 
emission of a variable electron population given by 



N^{E^,x) ^ A{Ee,x)E~°' X e"-B„ 



(13) 
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4 EXPANSION LOSSES 

In this work we use a uniform conical jet model with a constant 
opening angle and bulk Lorentz factor. This model represents an 
overpressured jet which is expanding ballistically. An overpres- 

4 

sured jet with a relativistic equation of state (p oc pa) expands 
freely with a constant opening angle and its expansion is essen- 
tially ballistic after a few scale heights ( [Begelman et al.| ( |1984[ l). 
A uniform conical jet with a constant bulk Lorentz factor is neces- 
sarily ballistic. This is because its radial velocity profile is v{r) = 
vo *" , so the kinetic energy of plasma associated with the radial 
expansion is constant along the jet, also, since the conical jet oc- 
cupies a fixed volume in space in the lab frame it does no work 
expanding its environment. This means that the expansion does not 
result in adiabatic losses in this ballistic model and so there is no 
conversion of the internal energy of the plasma into bulk kinetic 
energy. 

We will investigate the effect of including adiabatic losses on 
ajet in Paper II. 



5 SYNCHROTRON EMISSION 

We use the expression for the power emitted by an electron whose 
velocity is at an angle ^ to a uniform magnetic field of strength 
B{x) (see for example [Pacholczyk| (| 1 970|)) 



PiE„e,x) 






(14) 



For high energy electrons, synchrotron radiation is beamed 
along the velocity vector of the electron. In our model electrons 
are isotropically distributed in pitch angle, so integrating over 9 we 
find 



5.1 Opacity 

To calculate the observed emission through the jet fluid we need to 
know the opacity of the plasma. We calculate the opacity using the 
one dimensional radiative transfer equation. Neglecting scattering 
in a homogeneous plasma we find (see|Hughes|(|199l])) 



L 






(1 



(20) 



where I^ is the emitted flux, jv is the emissivity per unit 
volume, ki. is the opacity and D is the length of a section of jet 
plasma along our line of sight. Using the Rayleigh-Jeans limit and 
associating the thermal temperature of electrons emitting at their 
critical frequency with the kinetic energy of those electrons via 



kbT = Ee 



we obtain 



h 



2e2u2 



J_y_ 



kv{v,x) 



Ji'C 

1 5 

2e2 !/2 



(21) 



To simplify the geometry of the problem we assume that we 
observe at an angle smaller than the opening angle of the jet ^observe 
in the lab frame. We consider the emission from sections of the jet 
of width dx, where the opacity is a function of x. 



j^{v,x) = ne{u,x)P{u,x), 

2A{Ee,x)eaTB^/3^{ei^)^ 



3uKV,x) 



dx. 



3TTR{x)^dxfiom'^c^ 
j^{iy,x) = jo{Ee,x)u 2 . 

So using this in our expression for opacity we find 



(22) 
(23) 
(24) 



PBinglc(Ee,x) = aTCB{xfP^^'^. 



(15) 



Considering beaming effects and rotation frequency we make 
the simplifying assumption that electrons emit all their energy at a 
critical frequency i^c (this is a reasonable assumption if we integrate 
over a large range of electron energies) which is given by 



I^c{x) = 



3f^eB{x) 
4:Trm„ 



(16) 



see |Longair| ( |I994l l. Rearranging the above equation we can 
define the useful relations 



Ee = {€{x)iyc{x)y 



e{x)^ 



3eB{x) ' 



(17) 



We want to find the synchrotron emission from a slab in the 
jet at a distance x with width dx. The total power emitted at a crit- 
ical frequency is simply the number of electrons which emit at that 
frequency multiplied by the power they emit. 



k^{i', x) — 



jo{Ee,x)c^iy 
2eO-5 



(25) 



Using Equation[20]we see that if fc^D is small (the material is 
optically thin) the observed luminosity simplifies to the volume of 
the plasma multiplied by the emissivity as we would expect. If we 
observe down the jet within the opening angle we will see through 
to different depths at different wavelengths due to the wavelength 
dependence of the opacity. We see through to the layer where the 
total column optical depth is approximately 1 and beyond this there 
is an exponentially decreasing contribution from further layers. The 
total optical depth Ttot is the sum of the optical depths from each 
section. 

In order to calculate the optical depth of a section with lab 
frame width da;', for a photon travelling at an angle ^observe to 
the jet axis in the lab frame, we need to Lorentz transform its 4- 
displacement from the lab to the fluid frame, since we have cal- 
culated the synchrotron opacity in the fluid frame. We require the 
distance the photon travels in the fluid frame dr, this is related to the 
distance travelled by the photon in the lab frame by a simple Lorentz 
transformation. The 4-displacement in the lab frame is given by 



Ptot(j^, a;) = 
Ptot{i',x) = 



-Psingle(l^, x)Ne{u, x)dx, 

crTB'^l3'^A(Ee,x)e{x)(e{x)iy)^dx 
Spomic-^ 



(18) 



(19) 



This gives our total emitted power as a function of frequency 
in the centre of momentum frame of the electrons. 



da; 



cos(0|3 



/cdt'\ 
X' =i dx' \ = I dx' 

\ dy' J \dx' tanre'obst 



(26) 



We Lorentz transform this to the plasma rest frame to find the 
distance dr 
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' Tbuikdx' 



X 



\ dx'tan(6'obscrvc) 



")) 



dr = 7buikda::' 



cos(6lob 



/3 



(27) 



(28) 



For a section of lab frame width da;' the optical depth is simply 
the fluid frame synchrotron opacity multiplied by the distance dr 



Ttot{v,x) 



J X 



2:)7bull 



COS(^(7obscrvc/ 



13] Ax, (29) 



where we have converted the integration over dr to Ax using 
Equation 28 and the relation Ax = 7buikda; (due to Lorentz con- 
traction). An individual segment of width da; will emit according 
to 



1 5 

2e2 1/ 2 
Pv{x, Ax, v) — ■kR{x)Ax — (1 



^-k^i.,x)R(x)y (3Q^ 



We use the equation above to calculate the electron energy 
losses via Equation |ll| since this is the expression for the total syn- 
chrotron power emitted by a segment of the jet. We will observe the 
emission from each segment weighted by the total optical depth, 
which is equivalent to solving the radiative transfer equation with 
j.{z) = 0. 



Pi, = 'y^P,^{x,dx,i 



-Ttotil^tX) 



(31) 



The total synchrotron emission in the fluid frame, P^, is the 
sum of the emission from all the sections in the jet multiplied by 
fraction of radiation absorbed along the jet from a section. 



6 INVERSE-COMPTON SCATTERING 

In this section we calculate the inverse-Compton emission from an 
arbitrary electron population and arbitrary isotropic seed photon 
population. Let us consider the most general case of an electron- 
photon collision. In the lab frame we have an electron with energy 
Ee travelling in a direction, we are free to define our a;-axis (not 
to be confused with the jet a;-axis) to be aligned with the electron's 
momentum by a simple rotation. In these lab coordinates we have 
an incoming photon with energy E^ with its momentum vector at 
an angle of 9 to the negative x direction. We wish to calculate the 
scattered photon energy and angle in the lab frame after the colli- 
sion. 

To simplify the calculation we Lorentz transform into the elec- 
tron's rest frame (where the Klein-Nishina cross section is well de- 
fined) rotate our coordinates to deal with a simple stationary head- 
on collision, then apply a further rotation to realign our coordinate 
system with our original lab a;-axis and finally perform the inverse 
Lorentz transformation to revert to our initial lab frame. This is il- 
lustrated schematically in Figure [2] and is calculated using standard 
four-vector Lorentz transforms. 




L.T 




Rotation 




scattering 



m^c^ E/ 



•«- 



Rotation 
by -4) 



I 




Figure 2. Schematic diagram of the inverse-Compton scattering process 
defining our choice of variables. 



6.1 Obtaining the Inverse-Compton Spectrum 

We have calculated the scattering angle and energy of a photon 
which is scattered by an electron. It remains to relate this to an emit- 
ted flux. We approach this problem by discretely summing over in- 
terval values of all the dynamic variables (E^, E-,, 9 and (j)2 and x) 
and assigning a weight to each of the scattered photons which cor- 
responds to the number of photons produced per second in a section 
of width da: in the fluid frame. We then sum over the contributions 
from each section of the jet. 



Powers / E^sc!Ltt-v^eight{Ee,E^,9,(l)2,x) 

J E^,Ej,0,4,2,T 

xAEeAE.yA9A(l)2Ax. (32) 

6.2 Klein-Nishina Cross Section 

To calculate the weight function we need the cross section for the 
electron-photon interaction. The Klein-Nishina cross section is the 
first order QED cross section for electron photon scattering (see for 
example [Weinberg] (| 1 995)), it is given by 



Aa 
An2 2 



^a|sr,2p(£;,02)'(P«,<^2)+P(O2)-' ' '"-' 



1+COS (j>2) 

(33) 



P{E'^,M 



1 + 



; (1 +COS(f>2) 
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AQ,2 



TT Sin(/)2 



(34) 



where ops is the fine structure constant. The formula relating 
the number of collisions between electrons of energy E^ and pho- 
tons of energy E^ is 



Collisions(s ) 



N,{E.,x).'^^^^^.n,{E-,,e,x) 
xc(l - l3cose)AQ.2AEAE^AnAx, (35) 



n^i,(x,Ax) = 



Pv{x,Ax, ly) 



2nhh'R{x)Ax' 
and for the optically thick case (k^R(x) > 1) 

c3(e"B^ - 1) 
We define the average photon number density n"" (v) 



(39) 



(40) 



dfi^Trlsinelde, 



(36) 



where c(l — /? cos(^)) is the relative velocity of the photon and 
electron (see for example [Blumenthal & GouldH1970| l). Since the 
integral of the weight function is simply the number of collisions 
occurring per second we find 



Weight{Ee,E^f,e,(l>2,x) = Ne{Ee,X 



dcr(B;,( 



. nj{Ej,6, 
47r 



.c(l 



l3cos9)TT\sm(l>2\''r\sm9\. (37) 



We want to replace the integral in Equation [32] with a discrete 
sum over the variables. Under the assumption that the weighting 
function is approximately constant for sufficiently small parameter 
intervals in the sum, we can convert the integral over the weighting 
function into a discrete sum 

Power(W^) = ^ EjscM-V7eight{Ee,E^,9,(f)2,x) 

xAEedE.,AeA(j>2Ax. (38) 

This is our expression for the emitted power in the fluid frame. 
Using this method to calculate the emitted power has the advantage 
of calculating the inverse-Compton scattering exactly in the limit of 
small step sizes in the sum, whilst also obtaining the angular depen- 
dence of the spectral energy distribution. To do this we simply carry 
out the sum in equation[38]binning scattered photons depending on 
both their scattered angle, which can easily be related to their ob- 
servation angle, and their frequency. In this way we can calculate 
the emitted spectrum as a function of observation angle. 

In the next section we will calculate the photon number density 
n-y from the synchrotron emission. 



CM = E 



n^(z/, x)Ax 



(41) 



8 DOPPLER BOOSTING OF EMISSION 

So far we have calculated synchrotron and inverse-Compton emis- 
sion in the fluid frame. In our model the fluid is moving with a 
bulk Lorentz factor 7buik along the jet axis so the emission will be 
boosted depending on the Lorentz factor and the observation an- 
gle to the jet ^observe- We use primed variables to denote quantities 
defined in the lab frame and unprimed variables for quantities de- 
fined in the fluid frame. The lab frame jet opening angle O'^p^^-^^g is 
related to the fluid frame opening angle ^opening via 



7buik tan 6qpc 



= tan( 



(42) 



due to length contraction of the jet axis of the cone. The ob- 
served photon frequency in the lab frame is boosted by a Doppler 
factor (shown below) which is obtained by Lorentz transforming a 
photon from the jet fluid frame into the lab frame. The total ob- 
served power is related the emitted power in the fluid frame by 
multiplying by four Doppler factors one from Doppler shift in fre- 
quency, one from increased photon arrival times and two from trans- 
forming solid angle between frames. 



1 



7(1 -^COS0obsorvc)' 



(43) 



The observed power in solid angle dfl' at a Doppler boosted 
frequency u', P' [v'), can then be related to the synchrotron power 
emitted in the fluid frame at frequency v emitted into solid angle 
dfi, P{v), via 



7 SYNCHROTRON SEED PHOTONS 

At high observed photon energies(> IkeV), where inverse- 
Compton scattering starts to dominate the blazar spectrum, the jet 
is optically thin to photons. To calculate the inverse-Compton emis- 
sion from segments of the jet we need to calculate the photon num- 
ber density n^{x,v). When a section of the jet is optically thin 
(K^R{x) > 1) an infinitesimal section of width da; will receive ap- 
proximately the same amount of radiation from neighbouring sec- 
tions as it emits parallel to the a::-axis. We use the approximation that 
only the radiation emitted radially is unbalanced. So the volume of 
depth one light-second from the outer radial surface contains ap- 
proximately one second of emitted power. If the section is optically 
thick then the emission and absorption per unit length are equal and 
so the photon distribution in the object is that of a blackbody. Cal- 
culating n-y for the optically thin case we find 



P'{v') = P' 



= P{^) 



7(1 



P COS f/obscrvc^ 

1 



7(1 -/3cOs6lobserve) 



(44) 



(45) 



BL Lacertae has the very small cosmological redshift of 2 = 
0.0686, so we take the distance to BL Lacertae be d = 300Mpc 
and assume Euclidean space. 



9 RESULTS 

9.1 Modelling the SED of BL Lacertae 

In Figure [3k we show the results of our fluid jet model fitted by 
eye to the recent simultaneous muti-wavelength observations of BL 
Lacertae ( |Abdo et^T] pOTTJ). We stress how well our relatively 
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(a) A pure SSC fit to the simultaneous 2008-9 observations. 
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(b) This figure shows a comparison between models including and ne- (c) The electron energy distribution contained in a slab of fixed width one 
glecting electron energy losses light-second in the fiuid frame, shown after injection at the base and after 

losses at the end of the jet. 

Figure 3. Top: The observed spectrum of the model fitted to the recent simultaneous multi-wavelength observations of BL Lacertae. The model reproduces 
the observed spectrum at all frequencies Bottom: These figures show the effect of electron energy losses to the jet. The figure on the right shows that the high 
energy electrons have been depleted by energy losses at the end of the jet. The figure on the left shows that energy losses have the largest effect on the high 
frequency synchrotron and inverse-Compton emission of the jet, but do not affect the flat radio spectrum. 



simple model fits these observations: the model fits well to the data 
across almost all wavelengths using only inverse-Compton scat- 
tered synchrotron seed photons from the single population of jet 
electrons. There is a discrepancy at the lower end of the x-ray ener- 
gies, and we argue that this is likely from the accretion disc corona 
(see, e.g., |Jin et al.| ( |201 1) ). The optical/IR data is clearly not consis- 
tent with a smooth synchrotron component, but has been necessarily 
obtained with multiple instruments and may therefore suffer from 
absolute calibration offsets. By contrast, the optical-near IR spec- 
trum on BL Lac obtained by' Vermeulen et ar] ( |1995^ shows that the 
spectrum is clean and continuous in this region. 



Figures |3b and [3t show the effect of including electron en- 
ergy losses in the model. In Figure [3l we can see the evolution of 
the electron energy distribution along the jet. The figure shows the 
electron energy spectrum (contained in a slab of width one light- 
second) at the base of the jet immediately after injection and at the 
end of the jet after the slab has propagated through the entire jet and 
suffered maximal energy losses. We can see that losses are most sig- 
nificant for the highest energy electrons as we would expect, since 
the power of both synchrotron and inverse-Compton emission is ap- 
proximately proportional to 7^ of the electrons. The energy losses 
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Parameter 


Value 


Wj 


7.63 X lO^'^W 


L 


1 X IQiSm 


Bo 


3.63 X lO^^T 


Ro 


7.32 X 10"m 


A 

-^equi 


1 


^min 


5.11 MeV 


^max 


5.60 GeV 


a 


2 


e' 

opening 


9.7° 


^observe 


2° 


7bulk 


12 



Table 1. This table shows the values of the physical parameters used in the 
model used to fit to the BL Lacertae observations 



effectively lower the maximum energy cutoff of the electrons as 
they move along the jet. 

Energy losses are most severe close to the base of the jet due 
to its stronger magnetic field and smaller radius. This means that 
losses occur predominantly in a region close to the base which is 
much smaller than the total jet length. This results in the jet produc- 
ing a flat radio spectrum even when taking into account radiative 
energy losses since the electron distribution does not change signif- 
icantly along the vast majority of the jet. We find that the base of the 
jet dominates the synchrotron and inverse-Compton emission at the 
highest frequencies whilst the radio synchrotron emission, which is 
optically thick at the base, is dominated by the outer regions of the 
jet with larger radii and therefore a larger surface area for emission. 

In Figure [3b we can see the effect of the energy losses on the 
emitted spectrum of the jet. The figure shows that the low frequency 
synchrotron emission does not change significantly when losses are 
taken into account, however, the optically thin synchrotron emis- 
sion which is dominated by the highest energy electrons at the base 
of the jet is reduced. The effect of losses on the highest energy elec- 
trons also results in a reduction to the highest frequency inverse- 
Compton emission, for which they are responsible. 

From fitting to the observations of BL Lacertae we find a jet 
power of 7.63 x lO^^W, a maximum magnetic field Bq — 3.63 x 
10~^T at the base of the jet with base radius Ro = 7.32 x lO^^m 
and lab frame jet length L2 x lO^^'m. We also find a bulk Lorentz 
factor 7buik = 12, an electron energy index a = 2.0 and an obser- 
vation angle ^observe = 2°. We have fitted the observations with an 
equipartition jet ylcqui = 1 corresponding to energy being equally 
split between the electron population and the magnetic field. We 
have chosen the jet to be in equipartition because this is compatible 
with observations of radio lobes of Seyfert 1 galaxies being approx- 
imately in equipartition IMigliori et al.| ( [2007) ). 

In this model we have not included any reacceleration of the 
electrons through internal shocks in the jet. Internal shocks are a 
mechanism which can convert some of the large store of kinetic 
energy in the bulk motion of the jet into internal energy in the elec- 
tron population. This has been investigated by a number of authors 
([Marscher & Gear ( 1985), Kirk et al. ( 19981, Kaiser et al. (20001, 
[Spada et al., ( ,2001; , ,Sokolov et al., ( ,2004; and .Jamil et al., 1 ,2010; ) 
as a way of replenishing radiative and adiabatic energy losses of 
the electron population in jets. We expect that our estimates of jet 
opening angle and bulk Lorentz factor are larger than if we had in- 
cluded a form of in situ reacceleration, since a larger jet opening 
angle and Doppler factor reduce the effect of energy losses on the 
electron population. We intend to investigate the effect of different 
forms of reacceleration in this model in future work. 



The paper presenting the observations of BL Lacertae ( |Abdo| 
|et al.| pOTTJ) uses three models to fit to the most recent observa- 
tions. The first model from |Finke et al.|(2008^ is a single-zone time 
independent spherical SSC (synchrotron self-Compton) model. The 
fitted model uses a large Doppler factor of 26 and is over an or- 
der of magnitude from equipartition, it also fails to reproduce the 
low frequency radio emission. The second model based on work by 
|Ciprini| ( [2008[ l is a time-dependent single-zone SSC it uses two dif- 
ferent spherical blobs to fit to high and low frequency synchrotron 
emission and it fits well to the observations. The two blobs have 
Doppler factors 10 and 6.5 and the average equipartition value is 
Aequi = 2.8. The final model used is based on |B6ttcher & ChT] 
|ang| ( |200^ and is a single-zone time-dependent model with ERC 
(external radiation Compton) and also fits well to the data. The fit 
has a Doppler factor of 15, an equipartition fraction Acqui = 1.48 
but also includes an energy independent particle escape time of 
icsc = 60R/c, which is difficult to justify physically. Whilst these 
models are able to fit the high energy synchrotron and inverse- 
Compton emission they are not able to reproduce the observed flat 
radio slope of the spectrum, the self-absorbed spectra they produce 
are steeper. 

One of the major achievements of the conical jet model is to 
reproduce the radio emission observed at low frequencies. This can- 
not be reproduced in fixed radius, homogeneous, one-zone SSC, 
two-zone SSC and ERC models ( |Abdo et alT]p0I0| ) and |Marasc"hi] 
|& Tavecchio| l |2003) ) without introducing additional complications 
such as a two power law fit to the electron population or energy- 
independent diffusion of electrons. In fitting to the radio observa- 
tions as well as higher frequencies we gain additional constraints 
on our model parameters. In general we can gain the electron en- 
ergy spectral index from the slope of the inverse-Compton emission 
(see for example [Longair|(il994j ), the turnover of the flat spectrum 
at high frequencies then depends on the magnetic field strength 
Bo which is strongly constrained by the synchrotron turnover at 
~ IkeV and spectral index. Since blazars have such a variety of 
spectral behaviour with increasing variability at higher frequencies 
(see |Ulrich et al.| ( |1997[ )) it is surprising that most models neglect 
reproducing observations in the radio, since for most blazars the ra- 
dio is less variable and normally is observed to be close to a flat 
spectrum for a wide range of low frequencies. By not fitting to ra- 
dio observations, models reduce the number of constraints and lack 
a physical mechanism to reproduce one of the most stable features 
of blazars. 

Whilst our model has a number of free parameters we find 
in fitting the model to the observations that there are not large de- 
generacies between the different physical parameters. This is be- 
cause the different parameters affect specific aspects of the emission 
which are well constrained by the multi-wavelength observations. 
We find that increasing the jet power increases the emission across 
all frequencies without altering the shape of the spectrum signifi- 
cantly. The dependence of the inverse-Compton emission on both 
the synchrotron power emitted (SSC photons) and the total electron 
kinetic energy, which both depend on the jet power, leads to the 
inverse-Compton emission depending more sensitively on the jet 
power than the synchrotron emission. The electron spectral index 
a determines the slope of the optically thin synchrotron emission 
and the slope of the inverse-Compton emission which is reason- 
ably well constrained by the data. The magnetic field strength at 
the base of the jet Bo (which fixes Ro if Wj is specified) changes 
the frequency at which the synchrotron emission becomes optically 
thin and the power and peak frequency of the synchrotron emis- 
sion. The jet opening angle ^opening h^s a negligible effect on the 
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Figure 4. Top: The average number density of photons plotted against en- 
ergy for the model fitted to the simultaneous BL Lac observations. 
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high energy emission and only weakly effects the power emitted at 
a given radio frequency (smaller opening angles gives more radio 
power). The length of the jet L determines the lowest radio frequen- 
cies which are still optically thick to emission, this parameter is not 
well determined by the data since it has very little effect on the ob- 
servable emission. Finally, 7buik and 6'obscrvc together determine 
the Doppler boosting of the jet. The observation angle determines 
the optical path length travelled by synchrotron photons, but this has 
a relatively small effect, so it primarily changes the Doppler factor 
for a given value of 7buik- The bulk Lorentz factor of the jet deter- 
mines the Doppler factor but also has a more subtle impact on other 
jet parameters such as the jet opening angle in the fluid frame and 
the initial energy density in plasma in the fluid frame (see Equation 
[2](. This means that the bulk Lorentz factor has a significant effect 
on both the observed power and spectral shape of the emission. 

Figure H] shows the seed photon population averaged over the 
jet n"" (see equation 41 1 for the models in Figurelsl Fundamentally 
an ERC model with an unconstrained seed photon population can be 
used to fit almost any high frequency observation. So it is important 
to explicitly show the seed photon population used when fitting the 
model to the data to ensure clarity and to confirm that the photon 
population is realistic. 

9.2 Producing a Flat Radio Spectrum 

The observation of a near flat radio spectrum is a characteristic of 
most blazars ( lAbdo et al.| ( |2010| l). We find that our model produces 
a flat spectrum which is largely insensitive to the physical parame- 
ters of the jet. 

Figure [5] shows the overall flat spectrum produced by the 
model. The dashed curves show the observed emission from com- 
ponents at different distances along the jet. The figure illustrates 
the classic explanation of a flat spectrum; that different components 
have different turnover frequencies which add to give a flat spec- 
trum. We believe that this is the first time this component expla- 
nation of a flat spectrum has been calculated explicitly for a con- 
ical jet orientated close to our line of sight. The figure shows that 
components closer to the base of the jet have turnovers at the high- 
est frequencies and dominate the high frequency synchrotron emis- 
sion (due to stronger magnetic fields and smaller radii) before their 
observed low frequency emission is absorbed exponentially by the 
outer jet material. 

We have found that the production of a flat radio spectrum is 
largely insensitive to the parameters used in our model although 



Figure 5. This figure shows the observed luminosity of the BL Lacertae 
model looking at a small angle to the jet axis. The solid line is the total 
luminosity which is flat over a wide range in photon energy. The dashed 
lines show the observed contributions from thin slabs at different distances 
X along the jet which add to produce the flat total distribution. 



the magnitude and frequency range of the flat spectrum do show a 
dependence. The most sensitive parameters are the overall jet length 
which affects the extent of the flat spectrum to low frequencies and 
the value of Bq which affects the departure from the flat spectrum at 
high frequencies. This means that a flat radio spectrum is intrinsic 
to the uniform conical jet model which conserves magnetic energy 
and requires no fine-tuning. 

Figure |6] shows the dependence of the synchrotron emission 
on the functional form of B{x). We can see that if the magnetic en- 
ergy is conserved along the jet (B oc R{x)~^) then the synchrotron 
luminosity is flat over a wide range of frequencies. If the magnetic 
field evolves to conserve the magnetic flux parallel to the jet axis 
(B{x) oc R{x)^^) then we no longer observe a flat radio spectrum. 
Since a flat or nearly-flat radio flux is a well known characteristic 
of blazar-type objects, this is evidence for the conservation of mag- 
netic energy in simple fluid jets. This result was obtained by |Kaiser| 
P006| l who investigated a variety of fluid jet models and found in 
agreement with us that a model which conserved magnetic energy 
produced a flat spectrum whilst a model where magnetic flux is 
frozen did not. However, his calculations were limited to the case 
of observations perpendicular to the jet axis. 

Magnetic flux should be conserved along the jet both parallel 
and perpendicular to the jet axis. We have found that observations 
seem to indicate that the magnetic energy is conserved in a sec- 
tion of the jet as it expands and that magnetic flux perpendicular 
to the jet axis is conserved, whilst flux parallel to the jet axis in- 
creases along the jet. In order to conserve flux parallel to the jet 
axis this means that either the magnetic field is oriented randomly 
on small scales and so the net flux parallel to the jet is zero or the 
magnetic field in the jet is mainly toroidal and not poloidal. This 
is an interesting result since dominant toroidal magnetic fields in 
jets are predicted by MHD simulations and theory ( |Lyutikov et al.| 
( [2005[ >, |McKinney| J2006^ and [Komissarov et al.| ( |2007| )) and are 
supported by polarisation measurements and Faraday rotation maps 
( |0'Suflivan & Gabuzda| ( |2009l ) and |Kharb et aLlp059t ). 
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Figure 6. This figure shows the dependence of the synchrotron emission on 
the functional form of B{x). We can clearly see that if the magnetic energy 
is conserved (B{x) oc R~^) along the jet then the synchrotron luminosity 
is flat over a wide range of low frequencies. If the magnetic field evolves 
to conserve the flux through the surface x = constant then we no longer 
observe a flat radio spectrum {B{x) oc R~^). Since a flat or nearly-flat 
radio is a well known characteristic of blazar-type objects, this is evidence 
of the conservation of magnetic energy in simple fluid jets. 



jet which conserves magnetic flux parallel to the jet axis does not 
produce a flat radio spectrum. 

In this first paper in the series we have described a flexible 
model for a continuous jet which includes a detailed treatment of 
the synchrotron and inverse-Compton emission from the jet and ra- 
diative energy losses to the electrons. For the first time we have 
attempted a full-spectrum fit to a blazar using a conical jet model. 
We find that BL Lac can be described by a ballistic conical jet in 
equipartition with realistic physical parameters without a need for 
significant in situ acceleration or external radiation. In the next pa- 
per we will consider these latter effects of the model in the context 
of high-jet-power Compton dominant blazars. 
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10 CONCLUSIONS 

We have re-visited the uniform conical jet fluid jet model for AGN 
jets and for the first time used it to fit both the synchrotron and self- 
Compton components of the jet spectrum. In particular we have 
shown good agreement with these components in the recent multi- 
wavelength observations of the relatively low-jet-power source BL 
Lacertae from radio to 7-ray . Our model and conclusions may be 
summarized as: 

• We demonstrate that a fluid model of an axisymmetric jet with 
a fixed opening angle and conserved magnetic energy fits the si- 
multaneous observations of BL Lacertae presented by |Abdo et^ 
pOTT] !. 

• Our model conserves energy and particle flux along the jet, 
this allows us to estimate physical properties of the jet and central 
engine such as the total jet power and magnetic field strength from 
fitting the model to a blazar spectrum. 

• We evolve the electron population dynamically along the jet 
taking into account energy losses from synchrotron and inverse- 
Compton emission. We calculate the inverse-Compton emission by 
numerically integrating the exact expression and we explicitly show 
the synchrotron seed photon population. Unlike previous investi- 
gations, we integrate synchrotron opacity along the line of sight 
through the jet plasma taking into account the emission and opacity 
of each section of the jet. 

• We show that a conical jet model which conserves magnetic 
energy produces the characteristic flat radio spectrum observed in 
blazars in agreement with |Blandford & K6nigl| JI979[ > and |Kaiser| 
( (2006| >. 

• Under our conical jet assumptions, we find that for a source of 
this jet power (~ 10^^ W), that there is no requirement for signifi- 
cant re-acceleration within the jet to explain the observed spectrum. 

• We find that the production of a flat radio spectrum does not 
require fine-tuning of the model parameters. We show that a conical 
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